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Abstract 

We analyze the infinite range Ising spin-glass in a transverse field 
T below the critical temperature by a one step replica symmetry 
breaking theory (1S-RSB). The set of n replicas is divided in r 
blocks of m replicas each. We present results for different values 
of the block-size parameter m.The spin operators are represented 
by bilinear combinations of fermionic fields and we compare the 
results of two models: in the four state (4S)-model the diago- 
nal S? operator has two unphysical vanishing eigenvalues, that 
are suppressed by a restraint in the two states (2S)-model. In 
the static approximation we obtain qualitatively similar results 
for both models. They both exhibit a critical temperature T C (T) 
that decreases when T increases, until it reaches a quantum crit- 
ical point (QCP) at the same value of T c and they are both un- 
stable under replica symmetry breaking in the whole spin glass 
phase. Below the critical temperature we present results for the 
order parameters and free energy. 

PACS numbers: 64.60. Cn; 75.30.M; 75.10.N 
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1 Introduction 



The Ising model in a transverse field is widely studied for being the simplest 
system of interacting spins that presents a quantum critical point(QCP)[T]. 
We will not discuss here the extensive literature on results for several versions 
of the model, but we will concentrate instead in the transverse-field Ising 
spin-glass model in the ordered region. The experimental realizations of this 
model are LiHo x Y\_ x F± compounds[2]. 

When the transverse field T vanishes, our restraint 2S- model for two 
states reduces to the Sherrington- Kirkpatrick (SK)model of a spin glassjH] 
that in the ordered phase presents a landscape of many almost degenerate 
thermodynamic states separated by huge free energy barriers. The existence of 
these multiple states reflects itself in the breakdown of the replica symmetric 
solution of the saddle point equations below the ordering temperature with 
the onset of the Almeida-Thouless instability^ |3] and its replacement by 
Parisi's replica symmetry breaking solution[6j. The question is if such a 
picture remains true when quantum mechanical effects are included or if the 
quantum fluctuations will be strong enough to cause tunneling between these 
free energy barriers. 

The nature of the spin glass phase in the transverse field Ising spin glass 
model has been the subject of controversial results that seem to depend on 
the representation of the spin operators, because in the calculation of the 
quantum mechanical partition function special tools are needed to deal with 
the non-commuting operators entering the Hamiltonian. The method more 
currently used in the study of short-range [7] and infinite range (HI El UD1 
ITTl IT2] spin glasses in a transverse field is the Trotter- Suzuki formula [T3] . 
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that maps a system of quantum spins in d-dimensions to a classical system 
of spins in (d + l)-dimensions, and it is suited to perform numerical studies. 
In referencejlOj.it was found by using the static approximation that there 
is a small region in the spin glass phase where a replica symmetric (RS) 
solution is stable, while in reference [TT] it is predicted without the use of 
the static approximation that the RS solution is unstable in the whole spin 
glass phase. The instability of the RS solution in the ordered region was also 
obtained in reference |12j by analyzing a one-step replica-symmetry breaking 
(1S-RSB) solution, and using only partially the static Ansatz. 

Another way of dealing with the non-commutativity of quantum mechan- 
ical spin operators is to use Feynman's path integral formulations 
and to introduce time-ordering by means of an imaginary time < r < j3, 
where j3 is the inverse temperature. More recently|I7| new conflicting results 
were obtained by using the formalism in ref.^3] with the static approxima- 
tion and it was claimed that the RS solution is stable in most of the ordered 
region. 

A still different functional integral formulation consists in using Grass- 
mann variables to write a field theory with an effective action where the 
spin operators in the Hamiltonian are expressed as bilinear combinations of 
fermions |18| I19j. In a previous paper we used a replica symmetric(RS) the- 
ory with the static approximation and we obtained that the RS theory is 
unstable)!] in all the ordered phase JH]- The advantage of the fermionic for- 
mulation is that it has a natural application to problems in condensed matter 
theory, where the fermion operators represent electrons that also participate 
in other physical processes, like the Kondo effect [20J. 
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In the present paper we extend our previous results [TS] into the ordered 
phase with a 1S-RSB formulation. One problem with the spin representation 
is that the spin eigenstates at each site do not belong to one irreducible repre- 
sentation S z = ±2> but they are labeled instead by the fermionic occupation 
numbers n a = or 1 , giving two more spurious states with S z = 0. We call 
this the "four states" (4S) model, and despite the presence of these two un- 
wanted states the 4S-Ising spin glass model describes a spin glass transition 
with the same characteristics as the Sherrington - Kirkpatrick (SK) model 
[3| in a replica symmetric theory. A way to get rid of the spurious states 
was introduced before by Wiethege and Sherrington [2~T] for non-random in- 
teractions and it consists in fixing the occupation number + by means 
of an integral constraint at every site. We refer to this as the "two states" 
(2S)-Ising model. 

In sect. 2 we analyze the 4S-Ising and 2S-Ising spin glass models in a 
transverse field, within the static approximation in a one step replica sym- 
metric breaking (1S-RSB) theory^ . The static ansatz neglects time fluctu- 
ations and may be considered an approximation similar to mean field theory, 
that describes the singularities of the zero frequency mode^3]. Numerical 
Monte Carlo solutions of Bray and Moore's equations indicate that the static 
approximation reproduces the correct results at finite temperatures [22] • 

When T = the static approximation reproduces the results obtained by 
other methods, in particular for the 2S-Ising spin glass model we recover SK 
equations [3]. The results in both models are very similar; they both exhibit 
a critical spin glass temperature T C (T) that decreases when the strength T of 
the transverse field increases, until it reaches a quantum critical point (QCP) 
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at T c , T c (r c ) = 0. The value of r c = 2\/2J is the same for both models and 
the 4S-Ising and 2S-Ising models are identical close to the QCP, where T c « 
| (r c — T). We obtained for both models that the replica symmetric solution 
is unstable jl] in the whole spin glass phase, in agreement with previous 
results with the Trotter-Suzuki method^2]- in Parisis theoryjU], the first step 
of RSB is achieved by dividing n replicas in r blocks of m replicas each, and 
taking the limit n — > 0, r — > while keeping finite the block parameter < 
m < 1. The block parameter m is considered to be temperature dependent 
and it is determined by optimization of the free energy. In our work we release 
this constraint and analyze the results for fixed values of the block parameter 
m. We present results for the free energy, and order parameters that are 
consistent with ref . ^2] but not with that of ref . ^7] . Numerical solutions of 
the saddle point equations indicate that the RSB order parameter vanishes 
at T = 0. We present the model in sect. 2 and we leave sect. 3 for discussions. 
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2 The model and results 



The Ising spin glass in a transverse field is represented by the Hamiltonian 
H=-Y,Ji j S!S*-2TY l S? (1) 

ij i 

where the sum is over the N sites of a lattice and is a random coupling 
among all pairs of spins, with gaussian probability distribution: 

P(J tJ ) = e- J ^ N/wj2 ^N/lQiTjz . (2) 
The spin operators are represented by auxiliary fermions fields: 

s i = g fat -71 *!] 

Si = ^kW + a l°d ( 3 ) 
where the a\ a (a ic7 ) are creation (destruction) operators with fermion anticom- 
mutatiom rules and o =f orj, indicates the spin projections. The number 
operators n ia = a\ a a i(7 = Oorl, then Sf in Eq.® has two eigenvalues ±| cor- 
responding to riii = 1 — ritf, and two vanishing eigenvalues when = n^. 

We shall use the Lagrangian path integral formulation in terms of anti- 
commuting Grassmann fields described in a previous publication [TH], so we 
avoid giving repetitious details. We consider two models: the unrestrained, 
four states model and also the two states model of Wiethege and Sherrington 
where the number operators satisfy the restraint + = 1, what gives 
S- = ±|, at every site j2I] 

By using the integral representation for the Kronecker 5-function: 

5(n n + n n - 1) = i- jf* dxje^tor**-* (4) 
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we can express Z45 and Z 2 s, the partition function for the two models, in 
the compact functional integral form[T%] 

Z{fi = J D(<p*<p) J] i fj dx s e-*eW (5) 

where: 

rP d 

H(^ a (r)^ ja (r)} (6) 

and /ij = for the 4S-model while \ij = ixj for the 2S-model. We now follow 
standard procedures to get the configurational averaged free energy per site 
by using the replica formalism: 

1 Z(n)-1 

hm — (7) 



(3N n^o n 

where the configurational averaged, replicated, partition function < Z n > C)a 
Z(n) becomes, after averaging over J if 



a/3 



n d( i»p e 

-°° a/3 

IKTI^ r dx^e-^A^)} (8) 



with the replica index a = 1, 2, .., n, and 

Aj(q a p) = jD((p a (f a )exp[J2 A jar + 



pjHY,q«pSrSf] (9) 

a/3 
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We have introduced the spinors in the Fourier representation: 

/ 



(10) 



V 



and the Pauli matrices to write the spin glass part of the action in the static 
appr oximat ion , where 

and by introducing the inverse propagator 



we can write Aj a r 



7 1 = iuj + /ij + 



;i2) 



(13) 



We assume a one step replica symmetry breaking (1S-RSB) solution[U] 
of the saddle point equations by separating the n replicas in r groups of m 
replicas each : 



q + x, 



q a f3 = q + 5 
q a p = q 



\m J \m J ' 



(14) 



where q is the spin glass order parameter, x is related to the static suscep- 
tibility by x = Tx and 5 is the RSB parameter. The notation / (K) means 
integer part. 

The sums over a in the spin part of the action produce again quadratic 
terms that can be linearized by introducing new auxiliary fields, what makes 
the integration over Grassmann variables straightforward [T5] . In the limit 
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n — > 0, N — > oo we obtain the result for the free energy per site for the model 
with 2{p + 1) states, p = or 1 

I 1 POO POO 

(3F P = ((3Jf{[x+q} 2 -[6+q] 2 +m[5 2 +25q]} / Dz log / D Zl [2K p (q,x, 

Z Til J —oo J —oo 

(15) 

where we introduced the standard notation Dy = e 2 -JL and 

^ V27T 

/OO 
Dfcoshv'A (16) 
-00 

with 

h = jfiqzQ + Jy / 2(x - 5)f + Jv^^i, (17) 

a = (/?/i) 2 + (/?r) 2 (is) 

The saddle point equations for the order parameters are: 

/oo I poo 

- D *T5n*5?L D * i r- 11 *-* (19) 

/oo I /-oo roc Rh 

- C2 ° TK^?F { /-~ C2lXr /- ^ 7S smh ^ } (20) 

/oo I roo roo Rh 

~ B2 °E^? /_ { L B « 7T inh ^ } " * (21) 

where h(£, z) is given in Eq.(fT7j) and we called: 

h= r D({ m cosh ^ +pv s ^ } (22) 

J-OO A /\2 

We obtain for the de Almeida-Thouless eigenvalue [I] and entropy in both 
models: 

/oo 
-00 



, at _ -1 rVfl 7^2 / 1 P x 



. I-00 DziK™- 1 JZo sinh y^Affi 2 



} 2 (23) 



s 
k 



rn 
P 



-^J) 2 {[x + q} 2 -[S + q} 2 + m[6 2 + 26q}} + - / Dz \og Dz x K t 

Z Tfl J —oo J —oo 



The first terms of the Landau expansion of the free energy in powers of q 
and 5 gives: 

(3F P = (3F° - X -Dmq 2 - X -D{\ -m)[q + 5} 2 (25) 
where the coefficient D is: 

D = {(3Jf{l-2{(3J)\^f] (26) 

and the integrals are to be taken when q = 5 = 0. Here we consider < m < 1 
a free parameter and the equations will be solved for different values of m. 
We observe that when m = 1 we recover the results without RSB and q being 
the SG order parameter, while for m = we recover also the RS result, but 
with q + 5 as order parameter in agreement with ref. [U]. The critical line is 
given by the solution of D — 0. 

The numerical results for the critical temperature T C (T) were shown in 
ref . jTH] . For large values of the transverse field F the 2S-model and 4S-model 
are undistinguishable. When T = 0, the equations (|15|) for the free energy of 
the 2S-model (p = 0) reproduce Parisi's results for the SK model 0. For both 
models our equations coincide with ref. [12] except for eq.(|19p. because we find 
in our formalism that the static approximation is self-consistent. Finally, we 
comment on the de Almeida-Thouless instability. The numerical solution for 
X AT in Eq. (j23|) confirms that \ AT = on the critical line T C (T) and X AT < 
for T < T c (r), indicating the instability of the replica symmetric solution in 
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(24) 



the whole ordered phase for both models and contradicting recent results in 
ref.^Tj.We show in fig. 1 the free energy and entropy for the 2S model(p = 0) 
with RSB and m = 0.8, F = 0.5 J. In fig. 2 we compare the free energy for 
T = 0.5J in the 2S model with for different values of the block parameter m. 
In fig. 3 we present the results for the order parameters q and 5 for different 
values of m. In fig. 4 we present the results for the order parameters q and 
5 as a function of temperature for both models in the case of RSB, with 
r = 0.5J and m = 0.8. 
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3 Discussion 



We apply a 1S-RSB theory[6j to extend to low temperatures our previous 
study [IB] of two quantum Ising spin glass models in a transverse field by 
means of a path integral formalism where the spin operators are represented 
by bilinear combinations of fermionic fields, in the static approximation. The 
n replicas were separated in r groups of m replicas each and we analyze the 
results for fixed values of the block size parameter < m < 1. This differs 
from Parisi's theory[6j where m is considered to be temperature dependent 
and it is determined by optimization of the free energy. In the unrestricted 
four-states 4S-model the diagonal Sf -operator has two eigenvalues Sf = ±~ 
and two vanishing eigenvalues, while in the 2S-model the vanishing eigen- 
values are suppressed by means of an integral constraint [|2T] . Regarding the 
de Almeida-Thouless instability^, we obtained before [TH] that the replica 
symmetric solution is unstable in the whole spin glass phase for both mod- 
els, contradicting the results in ref. J7j but in agreement with the results in 
ref.|lljand in ref. obtained with the method of Trotter and Suzuki[T3]. 
The results we present here for the solutions of the saddle point equations, 
shown in fig. 3 and fig. 4, seem to indicate that the RSB order parameter S 
vanishes at T = for all values of the block parameter m. The free energy 
, shown in fig. 2, is maximized for smaller values of m and this is consistent 
with the results of ref.jB] where it is obtained m = at T = 0. However, if we 
introduce the value 5 = in our analytic expression in eq. ()24|) for the entropy 
we would recover the RS result, in discrepancy with the numerical results of 
ref.jH] at T = 0. We believe the origin of this discrepancy lies in keeping the 
block parameter m fixed. We notice in fig.3 that for m = 0.2 and T > 0.2J 
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the order parameters q and 6 cross, that is S > q. We may conjecture thar 
for m = this behaviour extends to all temperatures T < T c , in agreement 
with the numerical results |6J. 

To summarize, we study the ordered phase for the Ising spin glass model in 
the presence of a transverse field by using a fermionic representation of the 
spin operators. The RS theory is unstable for T < T c jTH] and we study a 
1S-RSB theory where we hold fixed the block-size parameter m . We derive 
saddle point equations for the order parameters that agree with previous 
results obtained with other methods Our theory has the advantage of 
computational simplicity and permit us to clarify the importance of m, but 
it does not describe RSB at T = for finite values of m > 0.2. We conjecture 
that RSB is restored at T = only for m = 0. 
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5 Figure Captions 

fig. 1 Free energy and entropy for the 2S-model with RSB, m = 0.8 ,r = 0.5 J 
as a function of temperature 

fig. 2 Free energy for the 2S-model with RSB and m = 0.2(full line), 
m = 0.4 (squares) and m = 0.6 (triangles), as a function of temperature for 

r = o.5J. 

fig. 3 Spin glass order parameter q and RSB order parameter S for T = 0.5 J in 
the 2S-model.m = 0.2 (full line),m = 0.4 (squares) and m = 0.6 (triangles). 

fig. 4 Order parameters for m = 0.8, T = 0.5 J in the 2S-model (dotted 
line) and 4S-model (full line) 
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